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Abstract. We study signatures in the Cosmic Microwave Background (CMB) 
induced by the presence of strong spatial curvature prior to the epoch of inflation 
which generated our present universe. If inflation does not last sufficiently long to 
drive the large-scale spatial curvature to zero, then presently observable scales may 
have left the horizon while spatial slices could not be approximated by a flat, Euclidean 
geometry. We compute corrections to the power spectrum and non-gaussianity of 
the CMB temperature anisotropy in this scenario. The power spectrum does not 
receive significant corrections and is a weak diagnostic of the presence of curvature 
in the initial conditions, unless its running can be determined with high accuracy. 
However, the bispectral non-gaussianity parameter /nl receives modifications on the 
largest observable scales. We estimate that the maximum signal would correspond to 
/nl ~ 0.3, which is out of reach for present-day microwave background experiments. 
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1. Introduction 

The precise data which led to our present era of observation-driven cosmology have 
largely been taken from measurements of the temperature and polarization anisotropics 
of the Cosmic Microwave Background [HE]- These anisotropics were apparently sourced 
by large-scale fluctuations of primordial origin, which passed inside the causal horizon 
while the universe was dominated by a hot plasma of relativistic matter and radiation. 
The existence of these long wavelength fluctuations can be inferred today by studying 
the imprint of the plasma oscillations which they seeded and which were left intact when 
radiation subsequently decoupled from matter. 

The origin of these superhorizon fluctuations is uncertain. With present-day 
observations, their properties are indistinguishable from a precisely gaussian random 
field. It is possible they were synthesized during an early epoch of inflation, in which 
case the detectability of non-gaussian features is known to be a precise discriminant 
between the simplest model, characterized by a single degree of freedom [3j and more 
complicated models, characterized by many degrees of freedom [3]. (For example, see 
Refs. [HJ, [6], [TJ, [HI El [TO] . ) In addition to its dependence on the dynamics of the inflationary 
phase, the non-gaussian signal may be sensitive to the composition of the universe before 
and after inflation. Subsequent to the end of the inflationary era, its evolution has been 
explored by several authors pj], [121 EE EU ESI [IB]- On the other hand, if inflation 
persisted for some number of e-folds prior to observable scales leaving the horizon then 
the impact of any pre-inflationary physics can be expected to be quite negligible [T7] . 

If inflation did not last long enough to erase all memory of earlier times, what 
information could be expected to survive in the statistics of large-scale fluctuations? 
Chen et al. [18] studied the impact of excitations above the Minkowski vacuu and 
found them to have a distinctive momentum dependence. Moreover, under certain 
assumptions the traces of these excitations could be made large. More detailed 
investigations were subsequently carried out in Refs. [201 IZE]. Statistical features of 
this sort could probe an enhancement of the correlations predicted by conventional 
inflationary models due to perturbative effects, but would necessarily have limited 
reach to detect a coherent modification of the background spatial geometry. In the 
conventional picture, where inflation is taken to have lasted for many e-folds prior 
to horizon exit of observable scales, the spatial geometry is taken to be flat to high 
precision. If this assumption is relaxed, however, then the background may be different. 
To study such a scenario one requires a new set of propagators and vertices, which 
are the elements out of which perturbation theory is constructed. Predictions obtained 
in this way could quite easily be rather different to conventional results. In the case 
of both perturbative or non-perturbative enhancements the effect diminishes as spatial 
slices inflate and isotropize, unless the modification is associated with new physics at 



| When promoted to de Sitter space, the choice of boundary conditions which reproduces the Minkowski 
vacuum on small scales is usually known as the Bunch-Davies state [19] . 
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an energy somewhat higher than the inflationary scalej§| 

According to present ideas the pre- inflationary phase could have been characterized 
by strong spatial disorder, of which the simplest manifestation would be non-vanishing 
spatial curvature. Assuming dark energy to be non-dynamical, observation presently 
requires —0.0178 < Q k < 0.0066 at 95% confidence, where fl k is evaluated today and the 
constraints are obtained by combining the WMAP 5-year dataset with baryon acoustic 
oscillation and supernova data |2j. As observations improve the allowed range of Q k 
can be expected to diminish, but if it remains consistent with zero it will presumably 
be impossible to rule out a small positive or negative value. One might be led to seek 
more information in the statistical properties of the large-scale perturbations, but it is 
already known that curvature is hard to detect in the power spectrum alone [22]. For 
that purpose it seems worthwhile to explore whether the situation is better for higher- 
order correlations. Thus, both to decide what can be learned about the pre-inflationary 
initial conditions, and to determine whether curvature at the level \Q k \ ~ 10~ 2 can be 
probed by non-gaussian statistics, one would like to obtain the inflationary bispectrum 
computed with non-flat spatial slices. In this paper we focus on the case of positive 
spatial curvature, corresponding to S 3 spatial slices, since in that case the observational 
limit is more favourable, as discussed above. 

How large an effect should we expect? The longest observable wavelength in the 
CMB contributes to its dipole, although one must be wary in extracting data at low 
multipoles because of scatter due to cosmic variance. Nevertheless, as a simple estimate, 
a calculation of the comoving distance to the CMB dipole (or £ = 2 mode; cf. Ref. [23] ) 
gives the diameter of the last scattering sphere as 

° dt -2 C dX (1) 

where we have specialized to S* 3 slices and, as above, the fli are evaluated today (at 
time t )- The decoupling time is denoted t^ c and e p is the scale factor. For the purposes 
of making an estimate, let us assume that the present-day spatial curvature takes the 
WMAP5+BAO+SN central value (at 68% confidence) tt k « -0.0050 [2J. The remaining 
cosmological parameters can be chosen equal to the best-fit concordance model where 
tt A = 0.721, fi r = 8x 10~ 5 and the CMB decouples at a redshift z dc = 1100. This gives 
d,£=2 ~ 0.46. On the other hand, in these units, the radius of curvature of the spatial 
slices is set by the scale r = 1. It follows that the CMB dipole receives contributions 
from wavelengths which are a significant fraction of the radius of curvature. Where 
this is the case, we do not expect to be able to neglect spatial curvature in calculating 
microwave background observables. 

In this paper we estimate the observational signatures of this scenario. We study 
the correction to the power spectrum from single-field, slow-roll inflation assuming that 

§ If there is a sensitivity to new ultra-violet physics, then quantum excitations can have correlations 
which are not well-described by the standard formulation. Such modifications are not inflated away, 
but instead become imprinted in each new mode which crosses the horizon [TSJ [20l [21] . 
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inflation began while f4 was non- negligible, and then consider the implications for the 
three-point function of the curvature perturbation. Our result agrees with the known 
flat-slicing result [3] in the small-scale limit k ^> 1 (where k is the wavenumber of the 
mode under consideration), in agreement with the naive expectation that curvature 
should play a less important role as inflation proceeds. On the other hand for small 
k the fluctuations can feel the curvature of the S 3 slices, leading to terms which are 
larger than the flat M 3 prediction by a factor of 1/e, where e = <p 2 /2p 2 is the slow-roll 
parameter which measures the deviation of the universe from exact de Sitter space. 
This enhancement is driven by quantum interference among the fluctuations at horizon- 
crossing, and is maximized close to the equilateral configuration where all momenta 
participating in the three-point function have an approximately equal magnitude]]]] 

The layout of this paper is as follows. In $2] we compute the action for a scalar 
field coupled to Einstein gravity in de Sitter space with S 3 spatial hypersurfaces. 
After including the effect of small perturbations in the metric, this controls how small 
fluctuations interact and become non-linear as the expansion of the universe draws them 
across the causal horizon. In order to obtain a description of the non-linearities which 
are observable in the CMB bispectrum, it is sufficient to expand the action to third 
order in the perturbations. In £J3] we study corrections to the two-point function, and 
show that curvature plays essentially no role in determining its structure at tree level. 
In §4] we give the third order action and in £J5] we calculate the three point function. 
We then discuss the non-gaussian parameter /nl imprinted in this scenario and show 
that we have the expected limit for modes of small wavelength. Finally, we give our 
conclusions in 

Throughout this paper, we work in units where the H = c = 1 and the reduced 
Planck mass is normalized to unity, giving Mp 2 = 87rG = 1, where G is Newton's 
gravitational constant. The unperturbed background metric is 

ds 2 = -dt 2 + e 2p 7ii dx l dx j , (2) 

and 7y is the metric on the unit three-sphere. This choice means that the constant 
curvature k is normalized to unity, and to avoid clutter we will omit writing k explicitly. 

2. Computation of the action 

To study the evolution of small fluctuations around the time of horizon exit, one must 
couple gravity to whichever fields are responsible for driving inflation. Since the energy 
density in these fields dominates the energy budget of the universe by assumption, 
small fluctuations necessarily induce perturbations in the background metric. The non- 
linearity of gravity endows these perturbations with self-interactions which cannot be 
neglected. 

|| Therefore, we do not expect that this enhancement by 1/e would provide an explanation for the 
observation of large non-gaussianity reported by Yadav & Wandelt [24] . 
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To simplify the subsequent calculation, it is very convenient to take the perturbed 
metric in the Arnowitt-Deser-Misner or so-called ADM form [25], where 

ds 2 = -N 2 dt 2 + h lj (dx i + N { dt) (dx j + N j dt). (3) 

In this equation N and N 1 are the lapse function and shift-vector, whose function is to 
guarantee reparametrization invariance of the metric. They are not propagating degrees 
of freedom but are instead removed by constraint equations coming from the Einstein 
action; and is the perturbed metric on spatial slices, 

h^e^H^ + T^). (4) 

The quantity ( is known as the curvature perturbation. If gravitational waves are present, 
represented here by a transverse and traceless first-order tensor F^, then they should be 
included as corrections to the unit three-sphere metric 7^. We have implicitly adopted 
a spatial gauge in order to write Eq. In the present paper we focus solely on 

the properties of scalar degrees of freedom, which can be summarized in terms of the 
curvature perturbation, and set r^- = 0. This is adequate for the purpose of computing 
the bispectrum of (. In the absence of perturbations, the lapse function and shift-vector 
reduce to unity and zero, respectively, and otherwise can be expressed algebraically in 
terms of the propagating degrees of freedom. 

The action is taken to be the standard Einstein-Hilbert term, together with a real 
scalar field 0, supplemented with the York-Gibbons-Hawking boundary term, 

S=^Jdtd 3 x v^jiV (R {3) -2V- h ij d^d^) +^ - E 2 + [</> - iV^] 2 ) j, (5) 

where Eij is a rescaled version of the extrinsic curvature, defined by = E^/N and 
given explicitly by 

Ey^^hij-ViNj-VjNi). (6) 

The scalar E is its trace, obeying E = E\, where spatial indices are raised and lowered 
with the metric h-ij. The covariant derivative compatible with this metric is denoted Vj 
and an overdot denotes a derivative with respect to cosmic time, so that x = dx/dt. 
There are two scalar degrees of freedom in this system: the fluctuations ( and 5(f). 
However, only one of these is a propagating mode; the other is a gauge degree of freedom 
and can be removed by a coordinate redefinition. Although a wide variety of gauge 
choices exist [26], two are of special interest. The first is comoving gauge, where the 
scalar perturbation is set to zero and only the curvature perturbation propagates, 

6(f) = 0, and h i:j = e 2p+2 S;- (7) 

Late on when a mode has passed outside the horizon 7^ can be effectively replaced by a 
flat metric on R 3 since the equations of motion are negligibly affected by the curvature 
of the S 3 slices. In this limit the scale of the metric is arbitrary and the equations of 
motion possess a symmetry under dilatations of the scale factor, where p t— > p + a for 
any constant a. It follows that on very large scales ( must be effectively conserved, 
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because it becomes indistinguishable from a shift in the background value of p. This 
conservation law has been shown to follow — independently of the theory of gravity in 
question [27] — from energy conservation and the absence of non-adiatbatic pressure. It 
is therefore convenient to express the answer in terms of ( when outside the horizon. 

Unfortunately, the calculation of correlation functions directly in comoving gauge is 
technically complicated. It is more convenient to work in the uniform curvature gauge, 
which is defined by 

5(fi = ip, and = e 2p ^ij. (8) 

In this gauge it is simple to compute the correlation functions of (p, but more work is 
needed to make a gauge transformation into the observationally-relevant variable £. It 
is this approach which we will follow in the remainder of this paper. 

2.1. Constraint equations 

The action (jHJ) does not depend on time derivatives of N or N\ and therefore its variation 
with respect to these fields does not produce evolution equations but only constraints. 
We need to solve these constraint equations to first ordej^j. The Hamiltonian constraint 
is 

-2V- jpiEyE* - E 2 ) - ^(0 - N^ 2 - h''o.,o, = 0, (9) 

where we have adopted the convention that a partial derivative with respect to x l is 
denoted with a comma, so that = dj<p. The momentum constraint is 

V, jl(£j - 8)E)\ - 1(0 - N%)^ = (10) 

In subsequent expressions, to simplify the notation, we will use <p = 4>{t) to denote the 
background value of the field and reserve <p for its perturbation. 

One can now specialize to the uniform curvature gauge. In this gauge, the solution 
to the constraints takes the form N = 1 + SN, where 5N obeys 

6N=\(£<p + x)\ (11) 



pV2 

and Ni = x,i, where \ is an auxiliary quantity determined by 

X = c(A + 3 - e)- 1 [1 (-L) + C J (12) 



.alt 

and A is the laplacian on the unit three-sphere. We note that both 5N and N l are 
suppressed by the square root of a slow-roll parameter, 

5N = 0(^) and N* = 0(y/e), (13) 

where e = <f) 2 /2p 2 is the slow-roll parameter defined in £QJ This scaling will be 
extremely important when we come to extract the terms of lowest order in the slow-roll 
approximation from the action. 

% As noted in Ref . [3] and given in more detail in Ref . [55] , the higher order terms in N and N l occur 
in the action multiplying the Hamiltonian and momentum constraints respectively, and hence they can 
be neglected. 
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2.2. The second order action 



In the uniform curvature gauge, Eq. (jHJ), the second order term in the action is 



S-2 



1 



dtd 3 x e 3 ^ [ - V V - - -y^ X + 2 - e- 2p Y J ^^j 

P P V 



-<p 2 + 0X<P + x 2 



(14) 



where x is given by Eq. ( TT2T) and A is again the laplacian on the unit three-sphere. 
Since Vq = —3e~ 2p + 0(e) and (|T3|) holds we have at lowest order in slow roll 



S 2 = ^J dtd 3 x e 3p ^ [0 2 + 3e~ 2 V - e"V'wl' 
The field equation for 99 which follows from this action is 

4 ( e 3p ^) - 3eV - e p A(p = 0. 
dt\ dt / 

As discussed below Eq. (JJJ), the terms depending on curvature scale out of Eq. (fl6l) in 
the limit e p — > 00 after which Eq. ffTBT) has a time-independent solution. 



(15) 



(16) 



3. Corrections to the power spectrum from curvature 

3.1. Scalar mode functions on the three-sphere 

We wish to solve Eq. ( fT6l) at lowest order in the slow-roll approximation. For this 
purpose it is first necessary to obtain an expression for e p at lowest order. One does 
this by Taylor expanding the scale factor about a particular time and collecting the 
terms of lowest order. In principle, the Taylor expansion can be performed at any point. 
However, because our final answer will involve an integral which receives its largest 
contribution from around the time of horizon crossing, we can achieve the greatest 
accuracy by expanding around the time of horizon exit. This is precisely analogous to 
the expansion of H and other background quantities around the time of horizon exit in 
the standard calculation on flat spatial slices. 

The background equations for the scale factor consist of the Friedmann constraint, 

3p 2 = \ft + V - 3e~ 2p , (17) 
together with the Raychaudhuri equation 

p = - l -^ + e- 2p . (18) 

In addition, the scalar field evolves according to the conventional Klein-Gordon equation 

<t> + 3p> + Vq' = 0, (19) 

which is not corrected in the presence of curvature. At lowest order in slow roll we can 
write the nth order time derivative of the scale factor in terms of alternating expressions 
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for odd and even n, 



d_ 

It 

d_ 

dt 



2n 



2n+l 



(20) 



(21) 



where an asterisk * denotes evaluation at the time of horizon crossing for the 
corresponding wavenumber k. It follows that after expansion around the time t = t#, 
the scale factor can be expressed at lowest order in the slow-roll approximation by 



e P ~ e P* cos h[(t - U)yjf% + e- 2 P*\ + 



P*e p * 



:sinh[(t - t*)y/ 'pi + e" 



•2p.l 



(22) 



J pi + 

We can see that this is the scale factor for a de Sitter spacetime with p and p matched 
to our slowly rolling background at t — i*. In fact, it is usually more convenient to 
work in terms of a conformal time variable, 77, defined by r\ = f dt' ' le p {t'\ where 

J Hhroat 

^throat corresponds to the throat of the de Sitter hyperboloid. In terms of the r\ variable, 
Eq. (1221) can be written 

1 



+ 0(e). 



A cos 77 

It follows that in the neighbourhood of r] 

p = p\Xsinr] + 0(e), 

where A = + t~ 2p * I pi and r] e (-7r/2,7r/2) for t G 
conformal and physical time is 



(23) 

77* the Hubble rate obeys 

(24) 

The relation between 



77 1 
tan — = tanh - 



(t - t„)/3*A + tanh 1 ( ^ 



2 2 L " ' - - VA/J (25) 

Combining Eqs. ( |T6l) and ( 1231) . we find that the equation for the modes at lowest order 
in slow roll is 

ipl + 2 tan 7] ip' k + (k 2 - 4 )</? fc = (26) 

where a prime ' denotes a derivative with respect to 77, and —(k 2 — 1), with fcGN- {0}, 
are the eigenvalues of the laplacian A. This is solved by the normalised positive 
mode [29] constructed to vanish at the south pole of the Hawking-Moss instanton [30] 
describing de Sitter space at horizon crossing, 

p\X(k 2 - 3) 1 / 4 I" ^ ( -1/2 + v/^3/2; -1/2 - y/W^3/2 

1/2 



a/2(A: 2 -4) 



sin 2 77 



{k 2 - 4) 
+ % „ sin 77 F 



3/2 



sin 2 77 



(27) 



3/2; -VW 
s/W^3~~"'~ V 3/2 

where F is the Gauss hypergeometric function and 7 2 = —1. The south pole is at 
77 = 700. This is equivalent to the usual formulation in which the positive mode is taken 
to approach its Minkowski value deep inside the horizon, where it is insensitive to the 
curvature of spacetime. 
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Ultimately we wish to use these mode functions to obtain analytic estimates of 
the n-point expectation values of fluctuations in tp around the time of horizon crossing. 
However, the presence of hypergeometric functions in Eq. ( 1271) means that it is awkward 
to use such expressions for this purpose. To obtain a more useful approximation, one 
may first re- write the equation of motion, Eq. (1261) . as 



Vk * _(A; 2 -5-2tan 2 77)^. (28) 



cos rj J cos rj 

The equation admits an "almost- WKB" solution, corresponding to 

^w-^L^costj, (29) 
y2k 



which becomes a good approximation for an increasingly large neighbourhood of rj = 
as k — > oo. Unfortunately, we cannot conclude that Eq. (1291) is sufficient to obtain 
an acceptable estimate of the three-point expectation value of scalar fluctuations near 
horizon crossing (where rj — > tt/2), because it freezes out to the wrong asymptotic value. 
This is a fatal deficiency. It would lead to an unreliable estimate of the magnitude of 
the bispectrum and therefore untrustworthy conclusions regarding the observational 
relevance of /nl- To find an approximate form for the mode which gives a good 
approximation over the entire relevant range of rj, we can use standard results relating 
to hypergeometric functions to find that at late times 77 — > 7r/2 the field freezes out to 

d _^ eWfc^/a ~ (30) 

(p _ 3) 1/4 ^ 2 (Jfe2 - 4) V2k* 

where the approximation holds for large k. It follows that if we take (pf to be given by 



^ 7k( C08r] ~k) e ' (31) 



then we obtain the correct value and derivative at rj = tt/2, and indeed this gives a good 
estimate for all rj. 

In Eq. ( 1311) and below, the label "cl" indicates that this set of mode functions form 
a good basis out of which we can build a quantum field when we come to quantize 
this theory. To do so, one introduces creation and annihilation operators at and ag 
which respectively create and destroy particles (as measured by an inertial observer 
deep inside the horizon) with momentum k. The quantum field corresponding to <p can 
be constructed by the usual canonical procedure, leading to 

f(v, *) = £ + Ql(x) <ff(v)a e ) (32) 

k 

where the S 3 harmonics are defined by AQ^ = —(k 2 — tyQt- The case k = 1 
is the homogeneous background and the case k = 2 gives purely gauge modesEl For 

+ To see that the k = 2 modes are pure gauge, suppose if2 satisfies Ay>2 = ~3(^2- It must also satisfy 
= — 7ijV2) as can be verified by checking each k — 2 mode, and therefore it is possible to perform 
a gauge transformation from the metric 

ds 2 = -N 2 dt 2 + e 2p 7 i3 (dx i + N i dt){dx j + N 3 dt), ip = tp 2 (33) 
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the purposes of studying fluctuations neither of these modes are relevant, and we will 
therefore always be interested in modes of wavenumber three or larger, so that k ^ 3. 

3.2. Appropriate coordinates 

We use coordinates in which the metric on S 3 takes the form 

ds\ = dx 2 + sin 2 X dfl 2 . (36) 

where dQ 2 , is the usual metric on S 2 . These coordinates are convenient because if we 
take our position to be given by x = then the CMB can be thought of as located on 
the copy of S 2 at x — Xl, where the subscript L denotes the time of last scattering. In 
this form the harmonics are given by 

Q s = n£(x)ir(M), (37) 

where k = (k,l,m), the Y™ are the usual normalised harmonics on S 2 and 



ni(x) = y ( t(/+ {2H - SW X^-'frra X)- (38) 

In Eq. (1381) . Cj^}^ is a Gegenbauer polynomial and k—1 / ^ |m|. The list k = (k, I, m) 
labels the quantum numbers associated with each harmonic, with k labelling the "radial" 
X-harmonic and (I, m) the usual harmonic labels on S 2 associated with the spherical 
harmonics Y™. Under the exchange m i— > — m, the harmonics have the propert}0 
that Q** = (— l) m <3_g. Therefore we may write the field (|32|) in the convenient form 

<pM = ^Qz(x)\vt(v)al + (-l) m vt*(v)aJ (39) 



3. 3. The two-point function 

From the above we arrive at the two-point function evaluated at late times 

In order to sum expectation values of <^ into expectation values of the curvature 
perturbation, (, one can use the non-linear 5N formalism, which will be described in 
more detail in §5.31 below. Here we merely quote the result for the two-point function in 

to a superficially different but gauge-equivalent metric 

ds 2 = -N 2 dP + e 2p 7y {d5? + fiTdfyidS? + N°di), (p = 0. (34) 

The coordinate transformation which allows us to pass between these two equivalent forms can be 
written 

t = t + fl, £ l = x l + AJI (35) 
* The (-l) m term here comes from [Y ; m (0, </>)]* = (-l) m Yf m (9, <p). 
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order to demonstrate that the power spectrum is consistent with the case of flat spatial 
slices, giving 

(feCfe) = (- 1 ) mi U > fe^a [i + °( k i 2 )} ■ (41) 

It follows that the power spectrum of fluctuations in ( is the same as in the R 3 case 
[22], up to very small corrections. 



3.4- The scalar spectral tilt 

In the absence of non-gaussianity, which will be studied in detail in §§1H5J the principal 
discriminant among models of inflation is the spectral tilt. Although Eq. (HIT) shows 
that the magnitude of the power spectrum of fluctuations approximately matches the 
flat space spectrum if k is not small, the tilt contains curvature terms from A and p\. 
These curvature terms must be taken into account when determining whether the model 
is compatible with the WMAP constraints on the spectral index n s . 

In the case of flat equal time slices one removes a term 1/k 3 from the two-point 
function when determining the tilt. This is because the tilt is a measure of the deviation 
from scale invariance. Since the degeneracy of modes on S 3 with a given k increases 
discretely with k, the notion of scale invariance is not natural on S 3 ; however, experiment 
proceeds on the supposition of flat equal time slices and thereby demands that we 
compute the tilt in the usual way. 

For large k the term p~ 2 e~ 2p * is small and so the tilt which follows from Eq. fj4T|) is 

n s - 1 = 277* - 6e* + 0(Ar 2 ) . (42) 
where the slow-roll parameters e and 77 are defined bjjjjl 

e = — and 77 = + — . 43 

2p z (j)p 2p 2 

The last term in Eq. fj42|) will be significantly smaller than the first two for k > 20, 
and so we will recover the usual result for flat spatial slices. In this case it is known 
observationally that the tilt is small [2j, so it follows that the combination of e* and 
77* in Eq. (142!) is too: generically, they are both small. More generally, Eq. (I42p shows 
that the two-point function is unlikely to be strongly sensitive to curvature unless good 
constraints can be obtained on a possible running. 

(j In the case of flat spatial slices, one often uses e H = —p/p 2 instead of the quantity e; here these 
differ by a curvature term, so that e = en + p~ 2 e~ 2p . The difference may be significant for those 
modes passing outside the horizon when there is substantial curvature, so we are taking Eq. (|43| to 
be the fundamental choice because the field is rolling slowly on a Hubble timescale. This also leads 
to (V') 2 /2V 2 w e/A 4 and V" /V w r\. Note that en need not be small for modes passing out of the 
horizon early on. 
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3.5. Projecting onto the sky 

In order to properly compare results in the S 3 case with those in the M 3 case one should 
project onto the sky and compare the resulting angular expectation values of (. For this 
we need to consider the S 3 harmonics evaluated at xl on the sky, where xl ^ 1- From 
Eq. (1371) and the eigenvalue equation AQ^ = — {k 2 — 1)Q% we see that H l k (x) obeys the 
defining equation 

-^( Sin2 ^ n ^) - ^^ n iW + (* 2 - mlix) = 0. (44) 
sin x "X ^ «X ' sin x 

In the limit x 1, A; ^> 1 and / 3> 1 the solutions become proportional to a spherical 

Bessel function, II^(x) ^ 3l(^x)- The correct normalization can be obtained by studying 

Eq. (138|) in the limit x and yields 

n*(x) « yf*jj(*x) («) 

for x "C 1, k 3> 1 and I > 1. In the limit xl ^ 1 the surface of last scattering is 
becoming close to our point of observation, compared with the radius of curvature of 
the S 3 . We therefore expect that curvature ceases to play any role and the two-point 
function projected onto the sky goes over to its flat space counterpart. Indeed, one can 
show that in this limit (and with the restriction I ^> 1 which guarantees we are looking 
on small angular scales) 

r°° dk 1 r i 2 

Jo k (piTrl J 

where the approximation of continuing the integral to is valid because the main 
contribution comes from the region k ~ Ixl ■ Eq- ( 146]) is equivalent to the well-known 
result in the R 3 case, which demonstrates the consistency of our calculation. 

4. The third-order action 

In the previous section we expanded the action to second order in the small fluctuation (p, 
which is sufficient to determine its two-point statistics and therefore the power spectrum. 
If we wish to go further, however, and determine the leading non-linearity then it is 
necessary to obtain a description of the process by which three (p quanta can interact. 
This information is provided by the third order term in the expansion of the action, 
which we will determine in the present section before going on to study the three-point 
function in £j5j 

After expanding S according to its definition, it follows that the third order term 
can be written 

S 3 = l - J dtd 3 x e 3 ^( - ^ V - V "5N V 2 - 5N^ 2 - 20iVV,* 

- SNe-^ipwj + 6p 2 5N 3 + ApA x 6N 2 - 5Nx\ijX lij 

+ 5N(A X ) 2 + 2j)8NN i <p }i - 4) 2 5N 3 + 2<p5N 2 <p) , (47) 
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where to avoid unnecessary clutter we have denoted the covariant derivative on the unit 
three-sphere by a vertical bar, so that Xu = VjX and = h^Xu. 

From this expression we require the leading order term in slow-roll. This will allow 
us to compute the expectation value (<p<p<p) to the first non-trivial order at horizon 
crossing, after which we must perform a gauge transformation to determine (CCO- After 
horizon crossing we expect the correlation functions of the curvature perturbation to be 
approximately conserved. 

To identify the leading slow-roll terms, we must relate the derivatives of the 
potential to the motion of the background field 0. For this purpose one can use the 
relations 

V™ =(j~ 3 ^) e ~ 2p + °( e3/2 )' and V o = ~ 3e ~ 2p + °( e )- ( 48 ) 

All other terms in Eq. f j4"Tl) are suppressed by at least one power of (fi/p, so the leading 
contribution to the action comes from V Q "'. On flat spatial slices this term is usually 
negligible [HI [32], since it is proportional to powers of (f> and <p. These terms are 
indeed present in Eq. (JHJ), but are accompanied by a term of order p/<p whose source 
is the curvature term in the background scalar field equation, Eq. (fT7|) . Accordingly, a 
significant three-point interaction can be present at early times when the scalar field is 
behaving in a way quite different from the flat slicing expectation. 

We are assuming that a slow-roll hierarchy exists at the time of evaluation, so it 
follows that the <fi contribution can be discarded and the leading contribution to the 
action can be written 

/3 
dtd 3 x^e p ^. (49) 


This is of order p/<fi ~ e -1 / 2 and is a qualitatively new contribution which is not present 
in the interactions among ip quanta on M 3 spatial slices [3]. This term is suppressed by 
two powers of the scale factor compared to the vacuum energy density and therefore 
appears in Eq. (|4"8l) proportional to e~ 2p , which implies that at late times where e p — > oo 
it no longer contributes to the interactions. Thus, for large k we can expect that the 
effect of primordial curvature disappears and we recover the standard flat space result. 

In £J5] we will determine the contribution that this interaction makes to the (£CC) 
expectation value. In fact, we will see that it gives rise to a term which behaves like 
e~ 2 , and which can therefore dominate over those terms which are common to flat 
hypersurfaces and S 3 hypersurfaces (to be described below), which behave at leading 
order like e -1 . Of course, as we have described, other terms in the slow-roll expansion of 
Vq" will dominate if e~ 2p is small — that is, if the curvature of the spatial hypersurfaces 
is small. We can only expect the term in equation (]4"9l) to dominate in the three point 
function if there has not been much inflation prior to the mode leaving the horizon. 

The next-to-leading order terms contribute proportional to e 1 / 2 . Their effect in the 
action can be written 

S 3 = I J dt d 3 x e^Vl^e-^p 3 + 3e~ 2p 5Np 2 - SNp 2 - 2pN l p\ t - 5Ne~ 2p p li p^ . (50) 
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These terms give rise to contributions which are common to both S 3 and M 3 , and are 
suppressed by a power of e compared to the leading curvature term (T49|) . They are the 
same as the terms obtained by Maldacena [3] . Adding both these contributions together, 
integrating by parts, and using the background equations for p and 0, we obtain 



and represents terms which vanish at late times. The term involving SL 2 /Sp is 
proportional to the leading-order equation of motion and therefore vanishes when we 
take the interaction picture field to be on-shell. However, it cannot simply be discarded; 
it records the contribution of boundary terms which were generated after integrating by 
parts and which we have not written explicitly [33]. The correct procedure is to make 
a field redefinition to remove these terms, by introducing a shifted field (p c , defined by 



This removes the nuisance terms proportional to the equation of motion and the 
boundary terms, giving a simplified action 

S 3 = - J ' dt d 3 x ^ (e? i p\ + e 5p <jxp 2 c (A + 3)" Vc) . (54) 

We may now compute the simpler correlation functions of ip c , and rewrite the result in 
terms of the interesting field p by using Eq. ( 1531) . 

5. The three-point function 

5.1. The ip correlation function 

The final step is to calculate the three-point function in the comoving gauge ([7j). As 
we have described, the comoving curvature perturbation ( is conserved after horizon 
exit in the absence of non-adiabatic pressure. Once a given scale falls back inside the 
horizon, ( can be used to seed the subsequent calculation of temperature and density 
fluctuations in the coupled baryon-photon plasma. 

To obtain the correlation functions of ( entails a number of steps. The first requires 
that we obtain the contribution from the reduced third order action (154"1) . and rewrite it 
in terms of the correlation functions of the original field (p using the field redefinition (152]) . 
Once this has been done it is necessary to determine how the correlation functions of 
( are related to those of ip. Fortunately, there is a simple prescription (the so-called 
"<5 N formalism") which is valid on large scales [3U [351 123 H]- The result is that ((t,x) 
evaluated at any time t later than the horizon-crossing time can be written 




(51) 



(52) 



<P = <fc + f(<Pc)- 



(53) 




(55) 
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where t* is taken to be a spatial slice on which the curvature perturbation £ vanishes, t 
labels a slice of uniform energy density, and N is the number of e-folds between these 
two slices. Note that this formula applies only in coordinate space and must be treated 
accordingly when transforming to harmonic modes. Eq. (|55|) can also be understood in 
terms of a further field redefinition which changes the action from the uniform curvature 
to comoving gauge j3]. The whole computation is performed with the interaction picture 
fields using either (l3Tj) or (1271) according to whether k is large or small. 

Let us first determine the contribution from the first term in the reduced third-order 
action, Eq. (|54l) . Provided we are only interested in tree-level amplitudes, the interaction 
Hamiltonian is simply given by minus the interaction term in the Lagrangian [36, 37J, 
so that ifi n t3 = — Ant 3- It follows that the contribution from the first term in (1511) . for 
large k, is 

/■7T/2-E 

(^k^k^k,) 2 -i dr) ([^(77/2)^(77/2)^(77/2), ^3(77)]) 

J 10a 

=■■■- dllr'^ 72 / £x ^ q ^ j+ c - c -' (56) 

where "c.c." denotes the complex conjugate of the preceding expression and k t = 
k\ + k 2 + k%. The function J is defined by 

j n /2 ~ £ sin 77 / i\( i\( i\ k 

J= / dr] — - — I cos f] — - — ) I cos rj — - — 1 1 cos r\ — - — ) e (57) 
J loo cos^V ki/\ k 2 ) V k 3 J 

In Eqs. ( 156|) - (|571) we have carried the integration over conformal time to within a small 
parameter e (not to be confused with the slow-roll parameter e) of future infinity. We 
can expect that this will be a good approximation because the integral receives its 
largest contribution from around the time of horizon exit and thereafter does not evolve 
appreciably, so there is little error in continuing the integration into the infinite future. 
Inside the horizon the integral is taken over a contour which turns a right-angle at 1] = 
and approaches infinity along the positive imaginary axis, corresponding to evaluation 
of the interaction in the Hartle-Hawking state. This is equivalent to the interacting 
vacuum of the full theory. 

Eqs. (|56|) - (|57|) give rise to a term in the (p correlation function which takes the form 

( ^^*3> 2 -^j|i7j / ^ ^Qk,Qk 2 Qk 3 ( - 2^ + ^-t^E M?) (58) 

In order to obtain the contribution this makes to the ( correlation funtion we will 
shortly see that the SN formalism tells us we must multiply by (— p*/^) 3 . The result 
will be proportional to e~ 2 k~ 8 p^, whereas the usual terms from the calculation with flat 
constant time hypersurfaces are proportional to eT x k~ % p\. So we expect the terms given 
in Eq. (|58|) to dominate for small k. 

On the other hand, some of the usual terms in the three point function are obtained 
by calculating the contribution from the second term in Eq. (l54j) . These terms are given 
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by 

Adding on the terms from the field redefinition (see section (15.21) ). given by Eqs. (1521) 
and fl53|) . and taking the dominant terms for large k, we obtain 

H*W 3 > 2 - / ^ VTQiQiQi {y E -^E fc3 + ^E (60) 

5.1?. Terms due to the field redefinition 

The field redefinition given by Eqs. (I52"l) - (l53l) contributes a number of terms to the three 
point function. These are calculated with the aid of a convolution (see Eq. (1631) ) and 
then contraction of the fields in such a way that only connected diagrams are produced. 
The first term of Eq. (1521) results in a calculation much like that which will be done in 
section §5.31 The second and third terms in Eq. (|52|) are given by a marginally more 
complicated calculation where one must integrate by parts to throw the (A + 3) _1 onto 
the Qj: in the convolution, and then contract the fields. The sum total of the redefinition, 
Eq. (El, is then 



(^ki ^kz ) c fci ^ c ls2 ^ c kz ) 

^ff^ / d* x SrQiQifii* (E^ 3 + 3 E^-E W) ( 61 ) 

5.3. The 5N formalism 

Finally, one uses the 5N prescription to obtain the full ( correlation function. 
Multiplying three copies of Eq. (|55p . taking correlations using Wick's theorem 
and truncating to tree-level terms which do not involve unconstrained momentum 
integrations, we find 

((^(kz) = (^) (n^k^ks) + \ {j£) [H^>*^)* 3 > + cychc] , (62) 

where * denotes a 'convolution,' 

V3 



{<p*<p) % (t) = J d 6 x^Qt(Z)^(t,x), (63) 

"cyclic" indicates that all cyclic permutations of {1, 2, 3} should be included in the sum, 
and N measures by how many e-folds the mode in question is outside the horizon (see 
description below Eq. ( 1551) ). Since we are working with a single-field model of inflation, 
the derivatives of N can be evaluated directly, yielding 
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5-4- The ( correlation function 



Collecting all these terms, the final result for the three point function evaluated on a 
late time slice is 



d 3 X y/lQ&Q&Qk 



~ 4 £^4£**4£^+#£** 



X 



i>j 



+ 



2pl 



12 



2 

i 



+ 6 ^-p£^ 2 -^:£ fc : 



(65) 



Of course, this uses the large k approximation for ip cl and so it can only be regarded as 
approximate. 

The leading term in this expression scales with momentum like [k~ 6 ]. This is the 
term computed by Maldacena, and is dominant for sufficiently large k, that is, on small 
scales. The correction terms we have computed scale with momentum like [k~ 8 ] . These 
terms can become dominant on larger scales. On sufficiently large scales, of course, 
one must remember that Eq. (1651) will be accompanied by other corrections that scale 
even faster as powers of momentum, like [fc~ 2n ] for any integer n ^ 3, and that these 
corrections will eventually overwhelm the ones we have computed. In this regime, for 
accuracy, one should use the full expression for the modes f[2T|) and perform a numerical 
assessment of the three point function. 



5.5. An estimate o//nl 

We are finally in a position to estimate the magnitude of the observable non-linearity 
parameter /nl which is produced by the sensitivity to curvature in Eq. fl65l) . We 
define the sign of /nl according to the WMAP convention, where the bispectrum is 
parameterized in harmonic space via 

B{h, k 2 , h) = ~/nl {P{h)P{k 2 ) + P{k x )P{h) + P(k 2 )P(k 3 )} . (66) 
5 

One must be careful in interpreting this formula for modes on S3, because when we 
work with the harmonics Q^, the fcj no longer have the same meaning as in flat space. 
The correct way to compare the magnitude of /nl between models with different spatial 
geometries is to project onto the sky and compare the non-linearity parameter in the 
angular expectation values of (. 

To obtain the correct projection, we begin by writing the three-C correlator in terms 
of a function k 2 , k 3 ) defined by 

(CfeCfcCfi.) = / d 3 x s/iQ^Q^iK h, h). (67) 

To obtain the angular correlation function of the £s, we follow the route outlined in §3.51 
the correlator in harmonic space is transformed back to coordinate space, and evaluated 
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at the radial distance corresponding to last scattering. In terms of the polar coordinate 
on S 3 , this is given by \ = xl- One finds 



where 



x ^m'xdx^xHlix^xmuk^h) (68) 
Jo 

crrr = / ^) c^^r^rc^ ( 69 ) 

and dfl 2 (x) is an element of solid angle on S 2 in the direction of the unit vector x. The 
integration over x is symmetric or antisymmetric about % = 7r/2, with the two regions 
interfering constructively for k t odd and destructively for k t even. Moreover, the H l k (x) 
under the integral can be replaced by spherical Bessel functions for the range of x-, k 
and I of interest. It follows that 

k ^ (A;iXi) } { 1 ^ 2 1 1 1 ^ 3 1 

x [l-(-l)^( kl ,k 2 ,k 3 ) (70) 

The influence of the [1 — (— l) fc '] term averages out to 1 as we sum over, say, k 3 and the 
summation over the ki may be replaced by integrals to give 

(Ci imi ( h m 2 ( lsms ) XL « C^T 2 2 7 Q 3 r 2 sin 2 X d X 

J 

dh k 2 j h (k lX )j h (kiXL)]{l <-> 2 Ul <-> 3k(fci,^ 2 ,fc 3 ) (71) 



'(ii+i) J I J I ■> 

First consider the part of £ which scales with momentum like [k~ 6 ] , coming from those 
terms in Eq. fl65|) which were computed by Maldacena and are present in the flat space 
expectation value. In order to carry out the integrations in Eq. ( 17TI) explicitly it is most 
convenient if the terms involving each ki in £ factorize, so that each ki integral can be 
evaluated independently. The necessary components of £ are 

^ kl)k2)h) = ^A — ( ^ yee + \y k * + \y ki e + 2 J^y k A + ...(72) 

There is a small variation with owing to the variation of p\, 0* and 0* with the epoch 
of horizon exit. Ignoring this slow variation, however, the only term which does not 
factorize involves l/k t . A similar obstruction was encountered by Smith & Zaldarriaga 
[38] . who found a factorizable form by introducing a Schwinger parameter, 

1 



dt e -t(k 1+ k 2 +k s ) _ ( ?3 ) 
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To understand which regions make significant contributions to the integral in 
Eq. (1711) one can account for the influence of the Bessel functions using a stationary 
phase method. In practice, one employs a WKB approximation to write these functions 
as an amplitude multiplied by a term with rapidly varying phase. An appropriate WKB 
solution can be constructed from the defining equation of the spherical Bessel functions, 

^ 2 [ Z m] = {^^-i)m Z )]. (74) 

The derivative of the WKB phase satisfies 



First consider contributions to the x integral in Eq. (JTTj) from the region where \ 3> Xl- 
There are two cases, depending whether k 2 is larger or smaller than 1(1 + l)/x 2 - For 
very large k, the stationary phase approximation implies that there is essentially no 
contribution to the integral. We are therefore left with the region k 2 < 1(1 + \)j\ 2 . 
Since x ^ Xl it follows that k 2 x\ 1(1 + 1), and therefore ji(kxL) is becoming 
exponentially suppressed. It follows that there is also a negligible contribution in this 
case. The only significant contribution comes from the region where x ^ Xl- In this 
region we can approximate sinx « x because Xl 1; the error this induces in the 
region x > Xl is immaterial because the integrand is exponentially suppressed there. 
Moreover, in view of what has been said about the region where x ^ Xl, there is very 
little error in extending the range of \ integration to infinity. Accordingly, we have 

, 2 \ 3 f°° 

(ChmiChmzChms/xL ~ C™ 1 ™ 2 ™* \~) / X dx 

J 

J ^dk x klj h (k lX )j h (kiXL)^ <- 2j|l ^ ^{h x ,h 2 >h)- (76) 

We are assuming that £ is factorizable, but if it is not a similar discussion applies after 
the introduction of Schwinger parameters. 

Let us relate this expression to the flat space limit, where xl <S 1. There is only an 
exponentially suppressed contribution to each k integral from the region ki G [0, li + 1], 
because if xl *C 1 then it must also be true that kfxt ^ h(h + 1) for each i, and 
it follows that the Bessel function ji^kiXi) is undergoing exponential suppression in 
this region. We therefore incur essentially no penalty in extending the lower limit of 
integration to 0, giving 

2\ 3 r°° 



(ChmiCl2m 2 Chm s )xL ~ ^T^T^i y~J \ X dx 







J dh klj h (k lX )j h (kiXL)}{l » 2j|l ^ 3 fc 2 , fcs) (77) 

This is the standard formula for the angular bispectrum of the curvature perturbation 
with flat spatial slices. It follows that in the approximate flat space limit, we can compare 
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£ between the cases of M 3 and S 3 spatial slices to find an estimate of the comparative 
magnitude of /nl- One finds 

5 r<£/4^.,., i Vfc3+ l Vfci2+ 2w 



NL 



+ 2 r — 12 



f4 



(7f 



The correction terms can be larger than the flat space result when k < e x / 2 . 



6. Conclusions 



We work in the scenario where the Universe is a slowly rolling but positively curved 
spacetime formed from a background with constant density spatial slices which are copies 
of S 3 , as opposed to R 3 . This possibility is compatible with observation if < 10™ 2 - 
10 -3 , depending whether the universe is taken to be positively or negatively curved. 
Indeed, interpreting present observational limits literally, the diameter of the CMB o^ =2 
satisfies approximately df=2 ~ 0.46 (in units where the radius of curvature is unity), so 
that the S 2 corresponding to the surface of last scattering spans about a twelth of the 
circumference of a great circle on the S 3 of constant time. On such a last scattering 
wall one can expect that the fluctuation eigenmodes on S 3 of k > 8 will contribute to 
the CMB. Modes of long wavelength, corresponding to small k, will pass outside the 
horizon early and therefore will be sensitive the presence of curvature. The question is 
to what degree a realistic observable such as the power spectrum, or / NL , can constrain 
the appearance of primordial curvature. We have shown that such modes, of small k, 
will not contribute significantly to the spectrum, but may contribute to /nl- 

We perform the appropriate selection of the vacuum for our slowly rolling de 
Sitter type spacetime, corresponding to the Hartle-Hawking state, and demonstrate how 
one should execute the calculation of both the two and three point functions on this 
manifold. Along the way we find exact expressions for the scalar fluctuation modes in 
global de Sitter coordinates, and since these are rather unwieldy we explore appropriate 
approximations. For a generic potential V one finds that the contributions to the action 
due to curvature at leading order in slow roll are remarkably simple. Further, we show 
how the S 3 slicing can be reliably compared with the conventional flat-slicing calculation 
on an S 2 surface corresponding to last scattering. In this way we can establish continuity 
with the M 3 calculation as one approaches the flat-space limit. 

We estimate the contribution of small- harmonics to the bispectrum. To compare 
with observations, it is most convenient to state the result in terms of Jnl, given 
in Eq. flTHj) . As written, this expression should be understood to be valid in the 
approximately equilateral case where all fc» have an approximately equal magnitude. 
However, we believe our calculation could be generalized, using methods similar to 
those of Maldacena [3], to accommodate the squeezed limit. 
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How large an /nl can be obtained? The k = 1 mode is the homogeneous 
background, and k = 2 is pure gauge. Therefore the interesting contributions to 
microwave background fluctuations must come from modes where k > 3. To estimate an 
upper limit, we specialize to the equilateral limit where ki = k for all i. After collecting 
terms in Eq. ( J7H1) . this yields 

5 145 
/N1 = --(23 £ ,-6„,)-^. (79) 

The first term, proportional to the slow-roll parameters e and rj, is precisely the flat- 
space contribution first derived by Maldacena. The second term is new, and accounts 
for the leading effect of curvature. It is irrelevant in the limit k —>■ oo, but can be large 
for small k. Since the signal is maximized by choosing k as small as possible, we can 
obtain an approximate upper bound by setting k = 3, leading to /nl ~ 0.3. It follows 
that the effect of curvature is marginally below the expected detection threshold for the 
Planck satellite, usually supposed to be of order / NL ~ 5 [39|, but might lie on the limit 
of what is practicable with futuristic technology such as a high-redshift 21cm survey 
jlQlST]. In the latter case, however, one would need some way to distinguish this small 
signal from the ubiquitous non-linearities of gravity itself [16] . 
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